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ABSTRACT 

To explain the observed neutrino masses through the seesaw mech- 
anism, a supersymmetric generalization of the Standard Model should 
include heavy right-handed neutrino supermultiplets. Then the neutrino 
Yukawa couplings can induce CP violation in the lepton sector. In this 
paper, we compute the contribution of these CP violating terms to lep- 
ton electric dipole moments. We introduce a new formalism that makes 
use of supersymmetry to expose the GIM cancellations. In the region of 
small tan /3, we find a different result from that given previously by Ellis, 
Hisano, Raidal, and Shimizu. We confirm the structure found by this 
group, but with a much smaller overall coefficient. In the region of large 
tan /3, we recompute the leading term that has been identified by Masina 
and confirm her result. We discuss the implications of these results for 
constraints on the Y v . 
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1 Introduction 

The discovery of neutrino mass has not only required revision of the Standard 
Model of particle physics but also of theories that go beyond the Standard Model. 
A compelling idea for the origin of the observed small neutrino masses is the see- 
saw mechanism. This requires the introduction of heavy singlet leptons, that is, 
right-handed neutrinos. In the context of supersymmetric theories, these singlet lep- 
tons belong to new chiral supermultiplets iVj, one for each fermion generation. The 
Yukawa couplings and soft supersymmetry breaking terms associated with these right- 
handed neutrino supermultiplets can play important roles in lepton flavor violating 
processes and in the production of the baryon number of the universe through 
leptogenesis 0. 

A particularly important aspect of this model is the appearance of new sources of 
CP violation. In addition to new CP violating parameters generic to new physics — in 
supersymmetry, for example, the phases of /i and the A — new phases are possible 
in the neutrino Yukawa couplings and in the neutrino B term (BMNN). Complex 
Yukawa couplings can lead to observable CP violation in neutrino oscillations, and all 
of these parameters can be the source of the CP violation that generated a fermion- 
antifermion asymmetry in the early universe |2l3j . 

To test whether the observed matter-antimatter asymmetry indeed arose from 
leptogenesis, it is necessary to determine the CP violating phases from microscopic 
measurements. There has been much analysis of CP violating observables in neu- 
trino mixing. In principle, it is possible to determine all seesaw parameters studying 
neutrino and sneutrino mass matrices but, in practice, it will be quite challenging, 
if possible at all, to extract all of these parameters in the near future [I]. Another 
possible experimental approach to test CP violation in the lepton sector is to measure 
the electric dipole moments (EDMs) of charged leptons j^]. There are in fact many 
possible ways that underlying CP violating couplings could give rise to lepton EDMs. 
Thus, it is important to classify these effects and, if possible, to learn how to separate 
them from one another. 

If CP violation is provided by phases of soft supersymmetry breaking parameters, 
it is straightforward to generate a contribution to lepton electric dipole moments in 
one-loop order. This possibility has been explored by many authors jH]. However, it 
is also possible to generate lepton EDMs in models in which the soft supersymmetry 
breaking terms are CP conserving, by making use of phases in the neutrino Yukawa 
couplings. A particularly simple context to study this effect is to consider models in 
which the soft supersymmetry breaking scalar masses are exactly flavor-universal and 
the A terms are exactly proportional to the Yukawa couplings. Such models arise in 
the simplest paradigms for gravity- mediated supersymmetry breaking [7j. The idea 
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of 'gaugino mediation' provides an attractive way to realize this scheme in the context 
of a complete unified or superstring model jH]. 

In this class of models with universal soft supersymmetry breaking interactions, 
the flavor and CP violation due to the neutrino Yukawa couplings is computed by 
integrating out the right-handed neutrino Ni superfields. This particular model has 
been studied in a number of papers, beginning with In particular, the contribu- 
tion of neutrino Yukawa couplings to lepton EDMs has been studied in this context 
by Romanino and Strumia JU] 5 Ellis, Hisano, Raidal, and Shimizu (EHRS) and 
Masina |T2] . These authors found that the analysis is complicated by GIM cancel- 
lations, so that the first nonzero contribution to the EDMs arises in two-loop order 
and has the form of a commutator of different combinations of the Yukawa matrices. 

The work of JT] and jT^j used a renormalization group equation (RGE) approach 
to evaluate the leading logarithmic contributions to the lepton EDMs. We thought 
that it might be valuable to extend these calculations by evaluating the complete 
contribution to the lepton EDMs without making leading-log approximations. In this 
paper, we present a new accounting method for the CP violating effects of the right- 
handed neutrino sector that makes this calculation straightforward. Our results, 
however, differ from those of [TT] even at the leading-log level. We confirm the 
general structure of the answers found by this group — in particular, the commutator 
structure noted in the previous paragraph. However, we claim that there are further 
cancellations not found in their papers that one must resolve to obtain the correct 
detailed formulae. We confirm Masina's result for the large tan f3 region, up to some 
minor factors, using a method that is much more transparent. 

The outline of this paper is then as follows: In Section 2, we specify the model 
in which we are working. In Section 3, we describe our procedure for integrating out 
the Ni supermultiplets and identifying CP violating contributions. In Section 4, we 
carry out this procedure for the leading CP violating contribution proportional to Y^, 
where Y u is the neutrino Yukawa coupling. We find a result that is parametrically 
smaller than that of EHRS by one power of a large logarithm. In Section 5, we 
reconsider the analysis of EHRS and show how that logarithm cancels out using their 
method. In Section 6, we give a formula for the lepton EDMs that arises from this 
contribution. 

In jT5j, Masina pointed out that, for large values of tan/3, a different contribution 
can dominate the evaluation of the lepton EDMs. This new term arises at one higher 
loop order, at order Y^Y^, where Yu is the charged lepton Yukawa coupling. In 
Section 7, we evaluate this contribution, which requires a nontrivial two-loop diagram 
calculation. 

In Section 8, we make numerical estimates of the electron EDM from our new 
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formulae and compare these to the results of other models of lepton CP violation. 

Our calculations in this paper look specifically at the terms resulting from inte- 
grating out the right-handed neutrino sector. We work from an initial assumption 
that the soft supersymmetry breaking terms are universal and flavor-independent. In 
a model with renormalizable interactions that violate the flavor and CP symmetries, 
this initial condition is not technically natural. Thus, there will in general be other 
CP violating contributions, for example, from the thresholds at the grand unification 
scale M GU t, that should be added to the formulae we present here. Because, in all of 
our formulae, the leading logarithmic behavior cancels due to a GIM cancellation, our 
terms are not parametrically enhanced over those from the GUT threshold. In spe- 
cific models, the GUT scale terms can be numerically smaller than the terms from the 
right-handed neutrino scale; the authors of j^j, for example, argue this for their SU(5) 
GUT model. In any event, our formulae are computed precisely for the effective the- 
ory of Section 2 with minimal subtraction (in the DR scheme) at Mqut- By noting 
this prescription, it should be straightforward to add GUT threshold corrections to 
our results when these are computed in a particular GUT model. 



2 The model 



We consider the super symmetric Standard Model coupled to three chiral super- 
multiplets iVj which contain the heavy right-handed neutrinos associated with the 
seesaw mechanism. The superpotential of the model contains the following terms 
involving lepton supermultiplets: 

W = Yi k e a pH la EiLjp - Y?e a pH 2a N t L^ - ^H la H w + ^M tJ N t N 3 . (1) 

In this equation, Ljp is the supermultiplet containing the left-handed lepton fields 
(ujL, £Jl)/3i Ei is the superfield whose left-handed fermion is £f L , and N{ is the su- 
perfield whose left-handed fermion is Vn. The iVj are singlets of SU(2) x £7(1). We 
introduce the right-handed neutrino masses My by hand, and we do not assume any 
a priori relation of these parameters to the other couplings in . 

Without loss of generality, we can choose the basis and phases of L, E, and N 
such that My and Y]? are real and diagonal. We will refer to the diagonal elements 
of these matrices as Mj, Y^. These choices exhaust the freedom to redefine fields, and 
so the matrix Y* 3 is in general off-diagonal and complex. The mass matrix of light 
neutrinos is given by 

ykiykj „ 



m. 
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If the neutrino Yukawa couplings Y* are of order 1, the requirement of small neutrino 
masses (m„ ~ 0.1 eV) leads to large values of the M k , of the order of 10 14 GeV. 

To the Lagrangian generated by (pj , we must add appropriate soft supersymmetry 
breaking interactions. In this paper, we will assume that slepton masses are universal 
at the messenger scale (of the order of Mgut) and that A terms are strictly propor- 
tional to the corresponding Yukawa couplings, with a real constant of proportionality. 
We will assume that the phases of fi and of the gaugino masses are zero. If these 
conditions are not met, it is possible to generate EDMs from one-loop diagrams, a 
possibility that has been exhaustively explored in the literature 

This restriction to universal, CP invariant, flavor invariant soft supersymmetry 
breaking terms is not a natural restriction of the model in the technical sense. It 
is violated by loop corrections due to the neutrino Yukawa couplings. In fact, our 
analysis in this paper is to calculate the CP violation induced by these corrections. 
Consequently, the effects we find can be cut-off dependent. As we have explained 
in the introduction, we will impose the universality and flavor symmetry of the soft 
supersymmetry breaking interactions as an initial condition, defined by minimal sub- 
traction in the DR scheme at Mgut- 

With this prescription, we will take the soft supersymmetry breaking terms for 
the lepton sector to be 

Lssb = — m l X! f*f~ m a^a^ a — \ Yiie a pHi a EiLip — Y* 3 e a p H2 a NiL j p) 
f 

-(^BvMiiNi) 2 + H.c.) - {^b H ^H x H 2 + H.c.) (3) 

where / collectively represents sfermions, and we assume that ao, bn and B v are real 
parameters. The parameters mo, m a , ao, and hn all have the dimensions of mass and 
are of order Msusy ~ 100 GeV — 1 TeV. CP violating phases arise both from the 
neutrino Yukawa couplings and from the neutrino A term, but, in this model, they are 
controlled by the same parameters. We should note that if any of the parameters ao, 
bu or B v has an imaginary part, the corresponding term can give a large contribution 
to lepton EDMs. This point is discussed in some detail elsewhere jSj. The specific 
effects of the B v term have been analyzed in 

In computing the effects of the Ni supermultiplets, it is convenient to work in 
components, keeping the auxiliary fields (the F fields) as independent fields. We use 
two-component notation for the fermion fields. With the effects of the Majorana mass 
term included, the propagators for the component fields of the Ni take the form 
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:^((?)n fc (-?)) = -^wf* ( 4 ) 

where a 11 = (1, <r) M and c = — icr 2 are 2x2 components of the Dirac matrices and 
the charge conjugation matrix. 



3 Radiative corrections due to Y v 

As it is well known, radiative corrections will distort the form of Eq. and 
break the exact mass degeneracy between the sfermions. In this section, we will 
focus on those radiative corrections to the parameters of Eq. Q that can induce 
CP-violating phase and EDMs, in particular, the effects of diagrams involving the 
neutrino Yukawa and A terms. We will discuss the form of the effective Lagrangian 
at scales just below the right-handed neutrino mass scale. When we compute the 
induced EDMs in Section 6 and 7, we will need to take into account some additional 
effects that come from renormalization group running down to the electroweak scale. 
In our analysis, we will always assume that the right-handed neutrino masses 
are much larger than the supersymmetry breaking mass terms, of order Msusy ; so 
that any contribution suppressed by M SVSY /M k can be neglected. In this limit, the 
calculation that integrates out the right-handed neutrino sector divides neatly into 
a part that corrects the supersymmetric Lagrangian and a part that corrects the 
supersymmetry breaking perturbations. 

First, we consider the radiative corrections to the supersymmetric part of the 
Lagrangian. We begin by noting that, to a good approximation, we can neglect dia- 
grams that include vertices from the supersymmetry breaking terms. Except for the 
fj, term, all coefficients in the supersymmetric Lagrangian are dimensionless, while all 
supersymmetry breaking terms have coefficients with mass parameters of order 1 TeV 
or smaller. Therefore, corrections to the dimensionless coefficients from the super- 
symmetry breaking terms are at most of order of Msusy/^cj completely negligible. 
Corrections to the /x term are at most of the order of fib a /Ml, again, a negligible 
correction. 

The radiative corrections within the supersymmetric theory are strongly restricted 
by the constraints of supersymmetry. All component fields within supermultiplet 
receive the same radiative corrections. By the non-renormalization theorem j!4j . the 
superpotential receives no corrections. The result of this theorem constrains only 
the leading term in a Taylor series in external momenta, but, since these diagrams 
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are evaluated at external momenta of order Msusy, terms that depend on external 
momenta are suppressed by powers of Msusy/^fc and can be ignored. Then the most 
general effective Lagrangian obtained by integrating out the multiplets will have 
the form 

L eff = J d 4 dLi(l + SZlY'Lj + J d^OEiil + 5Z E ) ij Ej + J d 2 6W + H.c. (5) 

Since the Lagrangian is real-valued, the matrices (SZl)^ and {5ZeY must be Hermi- 
tian to all orders in perturbation theory. Note that while (5Zl) %:> receives off-diagonal 
corrections at the one-loop level, (SZe) 1 -' receives off-diagonal elements only at the 
two-loop level because E does not have any flavor number violating coupling. 

To generate a lepton electric dipole moment, we require a flavor-diagonal matrix 
element of an electromagnetic form factor to have an imaginary part ^3] . However, 
the radiative corrections from the supersymmetric Lagrangian, treated to first order, 
will be proportional to the matrices 5Z L and 5Z E - Since the diagonal elements of a 
Hermitian matrix are real, none of these corrections, acting alone, can induce a lepton 
electric dipole moment. This is an important constraint, which we will continue to 
follow through our analysis. 

The soft supersymmetry breaking part of the Lagrangian receives corrections pro- 
portional to the supersymmetry breaking parameters. However, the form is still quite 
constrained. The most general effective Lagrangian has the form 

L S sb = -{ml + SmDijLlZj - (m 2 + Sm^ijE^Ej 

-{a Q Y B 5ij + 5A lj )e a pH la EiL jP + H.c. (6) 

Since Lssb is Hermitian, (5m 2 ^)ij and (Sm-^jij must be Hermitian matrices to all 
orders in perturbation theory. The A term can in general receive non-Hermitian con- 
tribution. However, we will show in Appendix A that, up to order Y v , , the corrections 
to the A term have the form 

5A ij = aoY^SZ'i , (7) 

where SZa is Hermitian. Here again, the form of the radiative corrections as Hermitian 
matrices limits their ability to contribute to electric dipole moments. 

To work with the effective Lagrangian written in (0) and ©, it is useful to bring 
the lepton and slepton fields into a canonical normalization by rescaling by (1 + 
SZ)^ 1 / 2 . Then the superpotential becomes 

W = -[(1 + 5Z E y^ 2 fY u [{l + SZLY^^^ALi (8) 
and the soft terms become 

LssBeft = -[(1 + 5Z L Y 1/2 (m 2 + 5m 2 L )(l + 5Z l Y 1/2 YUl, 
-[(1 + 6ZeY 1/2 (™ 2 + Sm\){l + SZeY^eJEj 
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-a [(l + 5Z E )-^ 2 Y t (l + 5Z A )(1 + bZ L )- x l 2 y ^H^l^ + H.c. (9) 



One more step is needed. To identify the mass basis for leptons, we need to 
re-diagonalize the lepton Yukawa coupling. Decompose the coefficient of (JEJ) into a 
product of a unitary matrix, a real positive diagonal matrix, and another unitary 
matrix: 

[(i + 5z E yv 2 Y £ (i + sz L )- x % = [(i + sv) T y £ (i + 5u)] {j (10) 

Then (1 + SV) T can be absorbed into the superfields E and (1 + SU) can be absorbed 
into the superfields L. The soft supersymmetry breaking terms now take a form 
similar to (|B|) : 

L S sb = -(ml + ArnfyijLlLj - (ml + AmfyijEjEj 

-a yu(Sij + AZ?l)e a pH la EiLip + H.c. (11) 

where 

(ml + Am?) = [(1 + 5U)(l + 5Z L y l ' 2 (ml + 5m\)(l + 6Z L )-^ 2 {1 + 5U)- 1 ] 
(m 2 + Am|) = [(1 + 5V)(l + 5Z E )- 1 / 2 (m 2 + 5m|)(l + 8Z E )- x /\\ + 5V)- 1 ] 

a y(i + az a ) = a y e (i + 5U)(i + 5z L f/ 2 (i + 6z A )(i + bz L y x i 2 (\ + 5U)- 1 . 

(12) 



At this point, the only signs of CP- violation from the neutrino Yukawa couplings 
occur in the coefficient functions listed in (|12|). It is still true that the first two 
coefficient functions are Hermitian matrices with real diagonal elements, and that the 
diagonal elements of the A term coefficient are real through two-loop order (order 
Yy\ For the mass matrices, this result is obvious. For the A term an additional 
slightly technical argument is needed, which we give in Appendix B. 

This implies that, through order Y"^ 4 , we cannot obtain a contribution to the lepton 
electric dipole moments from any individual term in (jll|) . However, we can obtain a 
matrix with an imaginary part by taking the product of two different matrices from 
fllljl . For example, 



a = Im 



AZ A Am\ .. (13) 



can have nonzero diagonal elements. Since both matrices are Hermitian, this quantity 
can be written more illustratively as 

c <=h( [ AZa ' Aml 0« • ^ 
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1+ ~ ~ I - 

l L E , L ^ 

Figure 1: The general form of the diagrams contributing to the EDM of a charged lepton 
£. The photon line should be attached at all possible positions in the diagram. 



Note that to compute Cf to order of Y^, it suffices to calculate AZa and Am| to the 
one-loop level. Through two-loop order, this is the only structure in the theory that 
can contribute to a lepton electric dipole moment. 

At three-loop order, products of Am| with the other matrices in can give 
additional contributions of a new structure. A specific CP-violating quantity that 
will be important to us is 

Di = Im ((Am|) T m£ Am|).. (15) 

This quantity also has a commutator structure, as we will see in Section 7. It is 
smaller than (fPfj) by a factor of Y^/4tt. Nevertheless, as we will see in Section 7, this 
term can give the dominant contribution to lepton electric dipole moments in models 
with large tan/?. To obtain the contribution from this structure of order Y^Y^, it 
suffices to calculate Am| to two-loop order and Am? to one-loop order. 

We can be somewhat more concrete about how the structures C, and Di arise from 
Feynman diagrams. Contributions to the lepton EDM's come from diagrams of the 
general form of Fig.^ in which a right-handed lepton and is converted to a left-handed 
lepton through a photon vertex diagram. A lepton line runs through the diagram, 
and the matrices (|12|) appear as insertions on this line. By the arguments just given, 
we need to consider contributions with two separate insertions. The product (fT%|) 
comes uniquely from diagrams of the form of Fig. |2Ja), with the photon inserted in 
all possible positions on the lepton line. The product (fT3j) comes from diagrams of 
the form of Fig. Efb). In the latter diagram, the left-right mixing contributes the 
factor of ni£. We will evaluate these diagrams in Sections 6 and 7, respectively. 
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Figure 2: Diagrams giving the dominant contribution to EDM of charged lepton t (a) for 
small tan (3, (b) for large tan (3. 
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Nk 




H2 

Figure 3: Diagram giving the field strength renormalization of the supermultiplet Lj. In 
this and the next few figures, we treat F components as independent fields; the F terms of 
iVfc multiplets have the propagators (jlj). 



4 One-loop corrections 

To estimate the lepton electric dipole moments at order Y„, we should next com- 
pute AZ A and Am|. According to the arguments of the previous section, only the 
leading-order contributions are needed. To this order 

Am| = 5m\ - m\bZ L AZ A = 5Z A (16) 

The factor 5Zl is most easily computed as the one-loop correction to the Fl field 
strength. There is only one diagram, shown in Fig. H2 its value is 

(5Zr) ij = (Y ki )*Y kj I d Pe (17) 

where pe is a Euclidean momentum after Wick rotation. 

The factor 5Z A arises from the diagram shown in Fig. 0] The vertex marked with 
a heavy dot is an A v vertex. The value of the diagram is 

WpWe + md ■ { } 

The tensor structure is exactly the same as in (fP7j) . This fact is used in Appendix B. 

The matrix 5m 2 L arises from the four diagrams shown in Fig. |3J The first diagram 
has two A u vertices; the other three have supersymmetry breaking mass insertions. 
It should be noted that there is a contribution in which in inserted into the F^ 
propagator, which results from the mixing of F^ with iV through the Majorana mass 
term. The final result is 

d A p E 



a Y ei (5Z A yi = -a Y ei (Y u k rY u k i J 



(27T) 



ml + al ml m 2 M% 

pUpI + Ml) + (p% + Mlf pUpI + MlY 



(19) 
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Figure 4: Diagram giving the one-loop radiative correction to the vertex Ag. The heavy dot 
is an A, vertex. 



which is quite similar to (J 17)) and (|18j) . except that some terms appear with two 
massive propagators. The small difference in structure between (II 9j) and the earlier 
equations will be significant. 

As we have explained in Section 2, we regularize these diagrams by dimensional 
regularization and minimal subtraction at the scale Mgtjt- This gives 



52% 



{5m\f 



(4tt) : 



1 (ykiyyV 



(An) 2 
i 2 

^ykiyykj 



logoff + 1 



2m 2 



(47T) 



log 



M 2 

GUT 

Ml 



(4tt 



^y ki^ *ykj 



M 2 

log^ + 1 



so that 



Aml ij 



(4tt) ; 



+ l] + a^[log- 



+ 1] 



.(20) 



(21) 



and AZa = 5Za- 



Now some significant simplifications appear. First, in evaluating (|14j) . we can drop 
any terms in Am- that are proportional to the tensor structure of 5Za- Thus, we 
can replace 



Am\ ij 



(4vr) 2V ^ ' ^ y • (22) 

Second, after making this simplification, we can drop any terms in 5Z % a that are 
proportional to the structure (Y^ l )*Y^. In particular, we can change Mqut inside 
the logarithm to any other value that is independent of k. We then find 



a 



(47T) 



[Yo.Yi] 



(23) 
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H 2 H 2 





Figure 5: Diagrams giving the one- loop corrections to the supersymmetry breaking L mass 
term. The heavy dot is an A v vertex; the marked insertion is a soft mass term uiq. 



where 

/If 2 

(y )« = (yo« = (y*W log j± . (24) 

As is explained just above, the expression for Cj actually does not depend on the 
parameter M^. It is convenient to choose to be the geometric mean of the 
to minimize the individual logarithms that appear in (|24j) . 

Our final result for Cj is simple and cutoff-independent. However, we remind the 
reader that this result is derived in the simple picture in which we ignore threshold 
effects at the GUT scale and regulate diagrams using the DR scheme. Because of 
the major cancellations that occurred in the simplification of C*, these threshold 
corrections, which depend in a model-dependent way on GUT-scale physics, can be 
of the same order of magnitude as (|23|). 



5 Comparison to the RGE approach 

It is remarkable that, to order Y^, the only contribution to the lepton EDM comes 
from the invariant C, and that there is no contribution from Im[A^]. This conflicts 
with previous results on lepton EDMs given by EHRS ^T] and Masina [T2j. In this 
section, we will compute the leading logarithmic contributions to Im[A^] using the 
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renormalization group method and demonstrate explicitly that they cancel. At the 
end of the section, we will compare our analysis to that of JT] and 



To carry out the renormalization group analysis, we must integrate the RGEs 
from an initial condition at Mqut to the heaviest N mass, M 3 , then from M 3 to M 2 , 
then from M 2 to M\. This procedure is valid only if M\ <C M 2 <C M 3 . Let us define 

t{Q) = j±p]ogQ . (25) 

and 

t 3 = i(M GUT ) - t(M 3 ) , t 2 = t(M 3 ) - t{M 2 ) , tt = t(M 2 ) - t{M 1 ) . (26) 

For a hierarchical spectrum of masses, we expect this procedure to reproduce the 
results of two-loop calculations up to the order Y^t 2 . 

It is very important to write the RGEs in such a way that the right-handed 
neutrino thresholds are accounted correctly. There are two aspects to this. First, one 
should, at each stage of integration, project out those iVVs that have masses above the 
scale at which the RGE is being evaluated. To discuss this, it is useful to introduce 
projectors P 3 — 1, P 2 — diag(l, 1, 0), P\ = diag(l, 0, 0), projecting onto the N 
mass eigenstates that are still active as we integrate through the various thresholds. 
Second, one should be careful to keep the matrices Yi and M diagonal, at least when 
heavy particles are integrated out. 

We found it surprising that it is necessary to worry about off-diagonal terms in 
M, and so we would like to illustrate this with an example. In the appendix of |12j . 
the RGE for the neutrino Yukawa coupling is given as 

AY 

^ = 3Y v Y^P a Y v + ■■■ (27) 

The contribution on the right-hand side arises from the diagrams shown in Fig. El 
The projector eliminates contributions from the right-handed neutrinos with mass 
Mfc > Q. Consider, in particular, integrating this equation down to a Q such that 
M 2 <Q < M 3 . Let t Q = t(M 3 ) - t(Q). Then the integration gives 

Y U (Q) = Y U (M GVT ) + 3Y v YjY v t 3 + 3Y u YjP 2 Y u t Q . (28) 

However, direct calculation of the diagrams in Fig. El with Euclidean external momenta 
with \Q 2 \ < Mf gives 

Y U (Q) = Y U (M GVT ) + 2Y u YjP 2 Y u (t 3 + t Q ) + Y v Y^Y„t 3 + Y u YjY u t Q , (29) 

since in the first diagram the contribution from A^3 in the internal line labelled N m 
has a propagator proportional to 1/{Q 2 + M 2 ) and so is suppressed for Q 2 <C M|. 
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Figure 6: Diagrams giving the terms proportional to Y 2 in the RGE evolution of the neutrino 
Yukawa coupling. 

The direct calculation is correct. The problem is that, at this level, the application 
of the renormalization group method is incomplete. The neutrino mass matrix also 
acquires off-diagonal terms from the RGE. For Q 3> M 3 

^ = 2(Y v Yj)M + 2M(Y„Yj) T + ■■■ . (30) 

Thus, we should, first, integrate all of the RGEs down to a scale of the order of M 3 , 
second, diagonalize the mass matrix M at this scale and rewrite the couplings in 
this new basis, third, integrate out its largest eigenvalue, and, finally, use the rotated 
couplings as the initial conditions for the stage of integration from M 3 to M 2 . 

It is not difficult to see that this prescription precisely eliminates the term that 
does not appear in (|29|) from (|28|) . The first diagram in Fig. |U] modifies the neutrino 
Yukawa coupling by a field strength renormalization factor Y v — > Z^^ 2 Y V . The equa- 
tion ()3U|) induces a similar modification in the mass matrix, M Z N 1/2 M(Z^ 1/2 ) T . 
The factors are the same in the two expressions due to the nonrenormalization 
theorem. When we now diagonalize M at the scale M 3 , the change of basis cancels 

— 1/2 

the off-diagonal 1-3 and 2-3 elements of Z N that affect Y v . 

To control this effect, one must integrate through all three thresholds by carefully 
solving the RGE for all couplings and mass terms. However, here we only wish to 
develop expressions for the effective couplings to order t 2 , in order to check the results 
of the previous sections. For this, it is easier and more direct to use the following 
procedure: First, we integrate the renormalization group equations for the couplings. 
Then we identify terms that correspond to diagrams such as the first one in Fig. H 
with decoupling internal lines, and we remove these contributions by hand. 

We should be careful also to remove diagrams with intermediate lines, since 
Fn also decouples, as we see from the last line of (@J). One- loop diagrams involving 
the supersymmetry breaking do term can produce mixing of L and Fl or N and F^, 
for example, as in the diagram shown in Fig. 0] In the RGE evolution of A u , we 
encounter a term in which the intermediate line is -F/v 3 . In this contribution, the 
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off-diagonal terms decouple and should be removed at the same time that we remove 
intermediate N 3 lines. Unless this is done, one cannot see the complete cancellation 
of Im[A«i] that we will present below. In an RGE analysis, this step requires treating 
Fn as a separate field in the Lagrangian and diagonalizing the N*Fn quadratic terms 
generated through RGE evolution. 

With this insight into how to treat self-energy terms in the RGEs, we can integrate 
the RGEs for coupling constants through the three thresholds. The renormalization 
group equations for coupling constants are as given the Appendix of [P2] . 

= W v K a + ■■■ 
= Y e K a + --- 
= AK a A v + hA v K a + ■■■ 
= 2Yt{Y v j*A v + A e K a + ■■■ 

= {m 2 L , K a } + 2{Y^P a m%P a Y„ + m\K a + A\P a A v ) + ■■■ 

= 2{m%Y}Y t + Y}Y t m\) + 4(y/m|F £ + m\Y}Y t + 4A e ) + ■ ■ • , (31) 

YjP a Y u and K a = Y u YjP a and the subscript a specifies the energy scale. 
Note that is a supersymmetry breaking mass and should not be mistaken for a 
large supersymmetric neutrino mass Mj. The terms not written explicitly in (j3*Tj) are 
terms with flavor structures such as Y v ■ tr\Y v Yj\ that do not contribute to EDMs. 
Using the prescription that we have explained above, we find different results from 
those found previously. When we solve for Ag and for Yg at a scale much smaller than 
Mi, we find that the imaginary parts of Yg and Agja^ are identical and are equal to 

Im [Yt{K 3 K 2 t 3 t 2 + KsKttzh + K 2 K x td x )\ . (32) 

Now one more step is needed. As in (J1(J|) . we need to choose a new basis for the 
leptons in which Y% is real diagonal after taking into account the radiative corrections 
due to the iV*. Since the imaginary parts of Yi and Ag/a found at the previous stage 
are identical, this completely removes the imaginary part of At, in agreement with 
our analysis in Section 3. 

It is not clear to us how the analyses of [TT] and [12] found nonzero diagonal 
terms of ImL4g] at this order. The discussion in [TTj does not discuss the issue of 
re-diagonalizing Yg and M. On the other hand, writes the renormalization group 
equations in a way that explicitly takes into account the decoupling of heavy states, 



dY v 
dt 

<m 

dt 
dA v 

dt 
dA e 

dt 
dm 2 f 

Lj 

dt 
dm% 

dt 

where K n = 
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as we have noted above. In addition, the calculations done in this paper keep Yi and 
M diagonal by adding terms to these RGEs following the 'rotating basis' prescription 
of Brax and Savoy [THj . The full RGEs considered are not written explicitly in |12j, 
but nevertheless they are used to generate the results that are quoted there [E| . One 
possible difficulty is that this analysis might not remove the terms with decoupling 
(N*F N ) propagators that we have discussed above (|3*T|) . 

The observation that lepton EDMs are proportional to the commutator of Y and 
Yi is the most important result of the analysis of EHRS [TT]. Once this result has 
been found, it is straightforward to obtain the correct order of magnitude for the 
contribution to lepton EDMs from the phases of neutrino Yukawa couplings. Thus, 
the qualitative dependence of EDMs on the underlying supersymmetry parameters is 
given correctly in this paper, even though the actual terms that produce the lepton 
EDMs are different. 



6 Electric dipole moments 



We are now ready to obtain the actual expression for the lepton EDMs by evalu- 
ating the class of diagrams shown in Fig. 



A general diagram of the form of Fig. ^ evaluates to the form 

1 



— iev (p')ca' J ' 1 ' 'q u u(p) 



(F 2i + iF, 



25iy 



(33) 



where F 2 is the usual magnetic moment form factor and i indexes the lepton flavor. 
The lepton EDM is then given by 



7 ^ S 
di — —eFo^i — 



1.9 x 10" u ecm) ■ F 25i ■ — ■ S 



(34) 



where 5* is the spin of the lepton and S = S/(h/2). 

We would like to find a contribution to the EDM proportional to Ci in (|T4*|) . For 
this, we should find a vertex diagram that depends on both Ai and m| and insert 
the flavor-violating corrections found in Section 3. The only such diagram is shown 
in Fig. EJa). The value of this diagram, as a contribution to F 2 + iF 25 , is 



iF 



2:> 



2^ 



Vr 



Ola 
1 



Vr 



Ola 



02(1 



dz j dx 



(A e — /i tan (3)m 2 u m a 
z(l-z) 2 



o 



(zm 2 . + (1 — z)(xm 2 g + (1 — x)m|)) 2 



(35) 
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In this expression, Vo is the unitary matrix that diagonalizes the neutralino mass 

b° = Y,VoiaX° a w° = Y,Vo2aX° a , (36) 

a a 

m a are the neutralino mass eigenvalues (with signs), c w = cos8 w , s w = sin8 w . 

The renormalization-group running of the soft supersymmetry breaking masses 
from the GUT scale to the electroweak scale corrections gives large but flavor-inde- 
pendent corrections to the E and L masses proportional to the GUT-scale gaugino 
masses. These terms do not contribute to the flavor-violating effects that give the 
dipole matrix element an imaginary part, but they should be taken into account in 
the denominator of (j33j) in evaluating this imaginary part. Thus, we have written 
f)35)l as depending on the electroweak-scale values of these masses m^ and m£. The 
full expression (I35|) can be checked against many papers on lepton dipole moments, 
for example, jl8llW2()j . 

Starting from (|33|). we replace A e by ci Yi5Za, and we include one mass insertion 
in the L line by acting on the integral with 

< 37 > 

A schematic version of this analysis for general flavor-violating perturbations is de- 
scribed, for example, in [21]. In our model, we take the indicated derivative of (JSSJ)? 
assemble the structure (AZ,tAm|), and replace the imaginary part of this object by 
(iCi) as given in ([T3)l. We thus obtain an expression for the lepton electric dipole 
moment of the form of (|34jh where 

2a ^ fV 01a \ fV 01a V 02 a\ mlm%a Q m a ([ Y , Yi])„ m\ m\ 

-T25i = T7-\E 1 1 16 : 9{—, —) > 

(47r) 5 V V c w J V c w s w J \m a \° % m l a m 2 a 

(38) 

where g(xL, xe) is given in Appendix C. For comparison with the results of the next 
section, we might write this result alternatively as 



25i 



(47r) 5 V V c w J V c w s w J \m a \ 6 ml ml 

M 2 

MiYfWiYryyr] ■ i-^i i og ^ . m 



There is a curious consequence of this result that follows from the fact that the 
trace of any commutator is zero. If this effect is the only source of the lepton EDM, 
we expect that 

d e /m e + d^/m^ + d T /m T = . (40) 
It is unclear to us how this simple formula could be tested to the required accuracy. 



18 



7 Electric dipole moments for large tan/3 



Masina ^2] has argued that, for large tan/3, a different contribution to the lepton 
EDM can be the dominant one. Looking back at the diagrams of Fig. Q and El studied 
in the previous section, we see that it is advantageous at large tan/3 to drop A e and 
keep instead the term /itan/3. We still need a second loop correction to combine with 
Am- , but this can come from inserting Am% in the right-handed slepton propagator. 

•1 

one of 

Yf m 2 T tan 3 /? 



Since Am|, arises at order YgY 2 , the new diagram has a size relative to the previous 



,„ tan/3 rs — — s— ■ , (41) 
(4vr) 2 P 8ttV sin 2 /? ' V ; 

where v = 246 GeV, assuming that the r lepton dominates the intermediate states 
in the matrix product. We will see in a moment that, whereas all large logarithms of 
Mqut cancelled out of the expression for lepton EDM in Section 6, the contribution of 
the large tan/3 region is enhanced by two powers of this large logarithm. As a result, 
the terms we will compute in this section can dominate over those we discussed in 
Section 6 for tan/3 > 10. 

To evaluate this contribution, we need to work out the mass insertion Am| in (jlljl . 
To begin, we must compute 5Z E and <5m|, up to the two-loop level. We need only- 
compute those two-loop diagrams that contain the maximum number of Y v vertices, 
since only these diagrams will contain factors of the CP violating phases needed for 
a contribution to the EDMs. 

Consider first SZ E . The one- and two-loop diagrams contributing to the field 
strength renormalization give 



5Z% = TYfA 



d A p E 



1 "^'J (2vr) 4 (p|) 2 

9Y (Y ki YY k 'jY ( dipE ( d ^ kE - (A9) 



(2k)* J (2KY( PE )[(k E - VE f + Ml}(k E ) 

The two-loop contribution of order YfY 2 comes from a diagram of the topology of 
Fig- Eta). Notice that the indices of 5Z E are transposed. This is appropriate because, 
in the figure, the direction of the arrows is reversed on the E lines. In addition 
to the two-loop diagram, there is a one-loop diagram involving the one-loop SZ E 
counterterm. This diagram has topology shown in FigE(b) and has the value 



where e = (4- d)/2. 
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Hi 

Figure 7: The structure of the two-loop contributions to Ani| that involve the neutrino 
Yukawa couplings. All possible particles from each supermultiplet should be put on each 
line of each diagrams: (a) proper two-loop contributions; (b) one-loop diagrams containing 
one-loop counterterms for Am- . 
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The contributions to 5m 2 E from one-loop diagrams and from two-loop diagrams of 
the form of Fig. [7^ a) are given by 



+ 



d d p E 2ml + a o 
(2n) d W 2 

d d p E f d d k L 



2Y ii (Y v ki )*Y u k % J 



2iryj {2irypl[{k E -VE) 2 + Ml)\k%y 
5mQ + 4<2q — 



V (k E -p E y + Mi) 

-2 M m^ J g ^ ^« + 4ag) , (44) 

Again, we only consider corrections involving the Y v that will contribute to the 
EDMs. The first line of (J44|) gives the complete one-loop contribution. The second 
line gives the two- loop contribution proportional to YfY 2 . This contribution comes 
from diagrams of the topology of Fig. Ufa). To compute 5m 2 E , we put E on the 
external lines and insert mg into the propagators or Oq into the vertices in all possible 
ways. The final piece comes from diagrams of the topology of Fig. Efb) with the 
counterterms associated with the one-loop corrections to Zl, m|, and Za- 

Note the order of the indices on (5m| J! in these contributions; this reflects the 
reversed direction of arrows on the external lines in Fig. [7| Also, note that the 
integrals over k E contain superpartners L, H\ that have masses of the TeV scale 
rather than the right-handed neutrino scale. These integrals are potentially infrared 
divergent, and we will replace (k E ) — > (k E + M| USY ) to regulate this divergence. 

To compute the final mass insertion Am|, we must now make the redefinitions in 
(112}) . If we expand in the Yukawa couplings, we find 

Am| = (6m% - m 2 6Z E ) + [SV , [8m\ - m 2 Q 5Z E ) } +■■■ , (45) 

where SV is introduced in JTUJ). To give a nonzero diagonal element in ()15p. we 
must expand the quantities in the first term to order YfY 2 . In the second term, 
we will obtain a nonzero contribution to (|T3|) by taking the one-loop expressions for 
5m 2 g and 5Z E together with the one-loop expression for 5V that follows from the 
5Z L contribution to JTUJ). Note, while the flavor-independent gauge corrections to 
5m 2 ^ and 6Z E commute with SV, the first terms in Eqs. (|42|) and (J44|) . although 
flavor- conserving, do not commute with 5V. That is why we have dropped the gauge 
correction in Eqs. (}4*2*|) and ()44j) while we have kept the Y 2 terms. According to the 
above equation, 

5VY 2 - Y 2 SV = Y t {6Z L )*Y t . (46) 

If we recognize that the one-loop expressions for 5m 2 E and 5Z E are proportional to 
Y 2 , we can use this expression to evaluate the second term of ([45)1 . Inserting the 
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value of 8Zl given in (jl7j) and transposing the matrix, we find a contribution of the 
same structure YfY},Y v Yi that we have in the other contributions to (Am|) T . 



Our final result for Am 2 f is then 



Am/ = 2Y £i (Y„ ki )*Y v k % 

d«k E d« PE 1 f^^^ 2mlMl 



(2iry (2rcY (k% + M s 2 USY ) 2 p|(fe - Pe? + M|) V (A* - p E ? + Mf 
d 1 1.1,. 2 2\ 

W? We + AW ' (Sp£ K + °' 



(4^iS )( (4^ (l0g ^if +1)) « +<1 ^- (47) 



The two-loop integrals are standard forms that are evaluated, for example, in the 
Appendices of and • Using these results, we find for the off-diagonal elements 
of Am| 

Am 2 / = j^yM^yy^ 



,2 , /,„2n 1 i„ 2 iw GUT , i__ GUT i__ IV1 K 



(Qm 2 + Aa 2 ) 



-2 g M| + ° g Mf ° g Mf USY 
flog^ur W" 

/If 2 M 2 M 2 



- 2m ° log ^- (3m » +a » )log SS (log ^f ■ (48) 

This formula is the exact result to order YfY 2 with ultraviolet regularization by min- 
imal subtraction at Mqut- It does not assume that the right-handed neutrino masses 
are hierarchial. The dependence on Msusy, with terms of at most one logarithm, is 
consistent with renormalization group evolution from the heavy neutrino scale to the 
weak scale. The leading logarithmic terms in this expression are in precise agreement 
with the result of Masina [T2*] . 

The dominant contribution to the EDMs for large tan /3 is now found by inserting 
both (Am|) T and Am| into the vertex diagram as shown in Fig.^b). The imaginary 
part of the diagram is proportional to the structure (|15|). The contributions to this 
formula have up to three powers of logarithms. Just as in the evaluation of Cj, we can 
take advantage of the fact that we are computing the imaginary part of the product 
of Hermitian matrices, which picks out the antisymmetric product of these matrices. 
In this case, the result contains the structure 

lm[(Y u k rY^ml(Yr)*Yr} , (49) 
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contracted by a function of Mk and M m . Note that the structure in is antisym- 
metric in the right-handed neutrino flavor indices k and m. When we antisymmetrize 
the expression contracted with this structure, the leading term with log 3 (Mq UT /M|) 
cancels out. However, while for C, the next subleading logarithm also cancels out, 
here it does not and so, unlike the previous case, the final result will depend on Mgut- 
More precisely, we find 

* (9 log log Mhz hg Mk + 9 log2 Mk log Mk 

ol yiog ^ log ^ log M| +Jlog M| log M ^ 

^rriT, M 2 , , M% , iWlx 

+ 6 log -f 2 1 08 W + 3 108 w + < 7 - * > log w) 

(9 log ^ log ^ log ^ + 9 log ^>l°g^ 

ML T] M 2 , , ., Ml Ml , M 2 , . .1/^ 

+14 loe — log — — + 51oe — - + 41oe — -loe — + (7 — 3ir 2 ) loe — - ) 

+ g M% g Ml + g Ml + g M 2 g M 2 USY + 1 J g M 2 J 



(50) 

The parameter is a mean right-handed neutrino mass. The precise definition of 
this mass is unimportant, because, as in (}2*3*j) . the various factors of Mn cancel out 
of ()50|) when we use the antisymmetry of the structure Im[Y^]. It is convenient to 
choose as the geometric mean of the Mk to minimize the individual logarithms 
in (f50j). 

If the right-handed neutrino masses are strongly hierarchial, as was assumed by 
|f 2| . (|50|) is enhanced by three large logarithmic factors. The formula we have given 
here is valid for any right-handed neutrino spectrum; for a spectrum without large 
hierarchies, the leading term still has two large logarithms. We also confirm the 
result of [Ej that the logarithmic dependence on Msusy cancels in the leading order 
of logarithms, though a small dependence does remain in a subleading term. The 
coefficient of our leading term is identical to that found by Masina. 



From this expression we obtain lepton EDMs of the form of ()34|) with 
8a fV 01a \ fV ia , V Q 2a\ ^m%m a tan/? 



25i 



(An) 7 V c w J V c w + s w J \m a \ 8 v 2 cos 2 /3 



23 




Hi <H!> 



Figure 8: Sub-dominant diagram contributing to EDMs. The vertices marked by heavy 
dots are Aj> vertices, and the marked insertion is a one- loop correction to m|. All of the 
momenta flowing through the indicated loops are of order MgusY- 



where h(xi,XE) is given in Appendix C. In the above formula, we have kept only 
the leading logarithmic terms, that is, terms with two large logarithms in the case 
of a general right-handed neutrino mass spectrum and with three large logarithms in 
the case of a hierarchial mass spectrum. If we wish to give an expression valid, in the 
general case, at the level of one large logarithm, we should include the corrections 
to Di from the TeV threshold, replacing the Msusy by the actual masses of L and 
H\. At the same time, we must include an additional contribution, shown in Fig. |HJ 
involving a two-loop integral with momenta at the TeV scale. An analysis to this 
accuracy is beyond the scope of this paper. 



8 Discussion 

In this paper, we have re-evaluated the contributions from neutrino Yukawa cou- 
plings to the lepton EDMs. In contrast to previous studies, we have shown that, in the 
mass basis of charged leptons, up to two-loop level, neutrino Yukawa couplings do not 
induce an imaginary part to the diagonal elements of the Ai term. However, complex 
neutrino Yukawa couplings can create EDMs for charged leptons through differences 
in the renormalization of the An terms and the slepton masses terms, through the 
diagrams shown in Figs. Q and |2l Our expressions for the lepton EDMs have the 



• (9m^ + 9a 2 m 2 + 24) (log 



^yY^mUY^yYr] h( 





(51) 
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same structure in terms of the neutrino Yukawa couplings as those previously given 
in jTl] and [T2*] . However, the form of the integrals contributing to F 25i is different 
because there is an extra mass insertion. Further, the overall size of the effect is 
decreased from the previous estimates, especially in the region of low tan (5. 

There is an important test for the origin of lepton EDMs in the neutrino sector. 
While complex a and /i induce EDMs both for the charged leptons and for the 
neutron, effects from the neutrino sector give zero EDM for the neutron while making 
nonzero contributions for the leptons. However, an EDM present only for leptons 
could in principle arise from an imaginary part to the Ai coefficient or the neutrino 
B term as well as from loop effects involving Y v . It is interesting to compare the 
magnitudes of the effects from loop level or tree level CP-violating contributions. 



In the low tan/3 region, the effect that we have computed in (|3a|) gives a lepton 
EDM of the order of magnitude 

9Q 4l M, 2 / 200 GeV\ 2 /m e \ . . 

4 -10-"^ log ^(^^ ) (-)ecm. (52) 

The effect from the large tan (3 region has a double logarithmic enhancement with 
respect to this value. If we estimate log 2 (Mq UT /M^) ~ 200, the effect that we have 
computed in ()51|) gives a lepton EDM of the order of magnitude 

di „ 10 -» (^f\ Y ; loM f^X) 2 {^] e cn, (53) 

These estimates can be compared to the current best limit on the electron EDM, 
d e < 1.6 X 10~ 27 e cm [21]. To achieve an EDM close to the current bound, we would 
need to have Y^\og{M\/ M^) ~ 100. However, the experimental limit on the rate of 
fi — > e7 places a limit on the Y u matrix elements |T] , 



M 2 

y.keyk, ^ ^ut < q j tan/3 , (54) 



so it seems unlikely to have such large values of the Y v . On the other hand, the effect 
of the neutrino B term leads to a potentially much larger estimate for electron EDM, 

in _ 27 Im(B„) ^ 2 ^ 200 GeV x 2 m u 

di ~ 10 — Y v — — e cm . (55) 

Msusy " V M SUSY J m e K > 

This could easily saturate the present bound. Also, we expect d^ ~ m^/m e d ei so if d e 
is close to its present bound, d M should also be observable in future muon storage ring 
experiments If complex Yukawa couplings are the only source of CP- violation 
and Y v ~ 1, the electron EDM should still be observed in the next generation of 
experiments, which aim for sensitivity to d e ~ 10~ 29 e cm 
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Over the longer term, CP-violating effects of complex neutrino Yukawa couplings 
can also be probed by lepton flavor oscillation in slepton production at colliders |2Tj . 
and perhaps also in sneutrino-antisneutrino oscillation [2H]. Better understanding 
of the systematics of leptogenesis can also play a role on constraining the neutrino 
Yukawa couplings. All of this information will complement the knowledge that we 
are gaining from neutrino oscillation experiments to help us build a complete picture 
of the neutrino flavor interactions. 
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A Corrections to the lepton A term 

In this Appendix, we will show that the contributions to the parameter SA 1 ^ in the 
effective Lagrangian (JHJ) up to order Yj are of the form of (J7J) with 5Za a Hermitian 
matrix. 

The structure (J7|) is easy to see. The A-term interaction must have one factor 
of Ye, with corrections coming from neutrino Yukawa couplings. Among the three 
fields coupled by Ye, only Lj couples to Y v . So the corrections we are must be 2-point 
diagrams attached to the Lj external leg. 

The difficult part is to argue that these diagrams add to a Hermitian matrix. 
For this point, we will give a formal argument and a diagrammatic argument. The 
formal argument is as follows: We can generate the universal A-term in our model by 
multiplying the superpotential JIJ by a new chiral superfield S, which we will eventual 
set equal to (1 + ao9 2 ). With this structure, the superpotential is not renormalized. 
However, a correction to the A-term of the structure (JJJ) can be generated by if loop 
corrections generate new terms in the Kahler potential of the form 

5K n = 8Z%(ES) n LiLj . (56) 

Then is a linear combination of the SZ n . However, since (|56|) is a correction to 
the Kahler potential, 5Z n must be Hermitian |2T?] . 
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Figure 9: A contribution to the renormalization of Ai in two-loop order from two one-loop 
diagrams. 



As an alternative to this argument, we provide an explicit diagrammatic anal- 
ysis up to order Y„. In particular, for the problem at hand, we are interested in 
polynomials in Y v that contribute to the Ai term and have a nonzero imaginary part. 

A given diagram with only one A" line could, in principle, contain structures Y v -Y*, 
Y v ■ Y u , or Y* ■ Y*. The vertex A? conserves the number of H 2 (in fact if 2 has no 
Ai coupling). However, the vertices A v and Y v change the number of H 2 by one 
unit. Since we ignore the masses of L and H 2 in diagrams involving N, the L and H 
numbers are conserved by internal propagators. Therefore, any radiative correction 
to Ai has equal numbers of Y v and Y*. 

Consider a diagram contributing to Ai with only one A" line and (n + m) Ye 
vertices. From the above result, we see that the most general polynomial that can 
appear in such diagrams is 

{Yj) n 52(y*)*Y?f(M k )(Yi) m (57) 
k 

whose diagonal elements are purely real. Notice that in the case of one-loop diagram 
shown in Fig. H] n — 1, m — and the matrix SZa [defined in flj|)] is Hermitian. 

Now let us focus on two-loop diagrams with more than one A" line. In such 
diagrams four A^-vertices are involved. A contribution from a product of two one- 
loop diagrams, as shown in Fig. El has the polynomial structure 

£ Y™ yj" fc * Y" k Y™ 13 * /1 ( M m ) / 2 ( M n ) . (58) 

mn 

It is non-trivial but easy to show that the functions fi and f% are of the same form. 
As a result, the matrix in is Hermitian. 

This brings us to irreducible two-loop diagrams contributing to A terms. The 
structures of these diagrams fall into three categories: 1) they can be of the form 

E YrY u mk *Y u nk Y^* 9l (M ni M m ) , (59) 
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2) they can be of the form 



Y^Yr(YryY u nk (Y u nk y 92 (M n , M m ) , (60) 

mn 

3) or they can be of the form 

£ YrY u mk (Y u nk nY u n rgs(M m ,M n ) , (61) 

m,n,k 

where gi, g%, g$ are real functions of M n and M m . The structure shown in (|HUJ) is 
manifestly Hermitian. If the functions gi(M m , M n ) and g^(M m ,M n ) are symmetric 
under M m «-> M n , the structures appearing in (jo^j) and (|dT|) will be Hermitian also. 
It is not very obvious that these functions have the required symmetry. But it is 
not difficult to show this by explicit examination of the diagrams. All the relevant 
diagrams are shown in Fig. Since the momenta propagating in the loops are of 
order of Mat, we can neglect the external momenta, which for our purposes are of 
order of Msusy- With this simplification, it can be seen that all these diagrams are 
symmetric under M m <-> M n . This completes the proof that, up to two-loop level, 
the diagonal elements of Ai remain real. 



B Expansion of A^l, eq. ( 1121 ). to order Y l 



In Section 3, we claimed that the diagonal matrix element 

1 + 8U)(l + 8Z L ) X ' 2 {1 + 8Z A ){1 + 6Z L )-V 2 (1 + 5U)- 1 ] 11 (62) 

is real through order . To order Y 2 , this is easy to see: The matrix element is the 
matrix element of the sum 

(8U + Uz L + 8Z A - Uz L - 6U) = 8Z A (63) 

and 5Z A is Hermitian. 

Working to order Y^, we first consider the separate contributions of order Y^ from 
each factor of 5Za, 5Zl, and 5U. The factors of SU cancel. The contributions from 
SZ A and 5Z^ are diagonal elements of Hermitian matrices and thus manifestly real. 

In addition, we must look at contributions in which two of these objects at a time 
are expanded to order Y 2 . To analyze these terms, we need the expressions for SZl 
and SZa given in (|2()j) . We also need an expression for SU. The definition of (1 + 8U) 
is that it diagonalizes the matrix 

(l + 8Z L )- l l 2 Y 2 (l + 8Z L )- 1 ' 2 . (64) 
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Figure 10: Irreducible two-loop diagrams contributing to Ag. The shaded boxes represent 
the full one-loop propagator corrections from the N, H2 supermultiplets. 
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Using first-order quantum- mechanical perturbation theory, we find that (6U)u = 
and, for i ^ j, 

Y 2 + Y 2 1 
1 h Y 2 -Y 2 2(4t:) 2[ v ] v 

Then any diagonal element of a product of any two of SZ L , 5Za, SU is of the form of 
the quantity 

{[Y^yy^og + 1]) ([Y?YYn\og ^ + 1]) (66) 
multiplied by a real- valued expression. No such term has an imaginary part. 



C Mass dependence of dipole matrix elements 

As we have explained in Sections 6 and 7, the dipole matrix elements that con- 
tribute to lepton EDMs contain derivatives of the expression 

J_ f{xL , = f 1 dz f 1 dx z(l-z) 2 

m\ ' Jo Jo (zm 2 + (1 — z)(xm 2 g + (1 — x)m|)) 2 

/■! z{\ - Z) I 1 1 \ 

Jo m 2 -, — m| \zm 2 + (1 — z)m 2 ^ zm 2 + (1 — z)m 2 E J 
where x E = m^/m 2 , xe = m E / m2 a - This expression evaluates to 

11 (I - x\ + 2x L \ogx L 1 - x\ + 2x E \ogx E \ 

f{XL - XE) = { — <T^y J ■ < 68 > 

To make one insertion of Am|, we need 

± g (x L ,x E ) = ^±f(x L ,x E ). (69) 

This has the value 

. , 1 /l - x\ + 2x L \ogx L 1 - x\ + 2x E \ogx E 

g(x L ,x E )- • 



2(x E -x L ) 2 \ (l-x L ) 3 (l-x E f 
1 ( 5-4x L -x 2 L + 2(l + 2x L )\ogx L \ 

+ 2(x E -x L )\ {i-x L y )■ 1 } 
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To make one further insertion of Am%, we need 

1 d 1 

—h(x L ,x E ) = — - — -g(x L ,x E ) . (71) 

This has the value 

,1 - x\ + 2x L logx L 1 - x\ + 2x E \o<gx E 

h(x L , x E ) - 



(x E -x L )z\ (1-xl) 3 (1-xb) 3 
1 /5 - 4x L - x| + 2(1 + 2x L ) logx L 



2(x E -x L ) 2 \ (l-x L ) 4 

5 - 4x g - x\ + 2(1 + 2x E ) \ogx E \ , 
+ (1-x^ ] - {72) 



For m|, m|, 3> we find 



1 



f(x L ,x E )tt- g(x L ,x E ) « -— = — h(x L: x E ) m , (73) 

2xl%e 2x l x e 6X l X e 

as we might have expected. 
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